
QUANTITATIVE INTERFEROMETRY IN THE SEVERE ACOUSTIC 
ENVIRONMENT OF RESONANT SUPERSONIC JETS 


Carolyn R. Mercer 
NASA Lewis Research Center 
Cleveland, OH 44135, USA 
E-mail: r.arolvn.r.mercer@lerc.nasa.goy 


Ganesh Raman 

NYMA Inc., NASA LeRC Group 
Brook Park, OH 44142 
E-mail: panesh raman@lerc.nasa.goy 


ABSTRACT 

Understanding fundamental fluidic dynamic and acoustic processes in high-speed jets requires 
quantitative velocity, density and temperature measurements. In this paper we demonstrate a new, 
robust Liquid Crystal Point Diffraction Interferometer (LCPDI) that includes phase stepping and can 
provide accurate data even in the presence of intense acoustic fields. This novel common path 
interferometer (LCPDI) was developed to overcome difficulties with the Mach Zehnder 
interferometer in vibratory environments and is applied here to the case of a supersonic shock- 
containing jet. The environmentally insensitive LCPDI that is easy to align and capable of measuring 
optical wavefronts with high accuracy is briefly described, then integrated line of sight density data 
from the LCPDI for two underexpanded jets are presented. 

L INTRODUCTION 

There is renewed interest in the study of supersonic jets due to advances in high speed jet 
propulsion, supersonic combustion, and jet noise suppression for the next generation supersonic 
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commercial transport. Understanding fundamental fluid dynamic and acoustic processes for these 
applications recjuires quantitative velocity, density and temperature measurements. Thermocouples 
and pressure transducers can measure these parameters at discrete points, but these devices are 
invasive and alter the flow that they measure. Optical methods such as schlieren flow visualization 
have been used to get real-time, qualitative data, but it is difficult to extract quantitative data from 
schlieren images. Laser velocimetry, particle image velocimetry and doppler global velocimetry are 
used for velocity measurements but require seeding. Rayleigh scattering can be used for planar 
density, temperature, velocity and species measurements. Interferometry also yields density data, but 
to date it has been difficult to make quantitative interferometric measurements because of the high 
vibration environments associated with supersonic jets. This severe acoustic environment includes the 
three components of supersonic jet noise (Fig. 1). One of them is a discrete tone referred to as jet 
screech (Powell 1 ). The other two are shock-associated noise and jet mixing noise, both broadband 
sources. The screech resonance involves growing hydrodynamic disturbances in the shear layer that 
interact with the shocks to produce acoustic waves. These waves then propagate upstream and excite 
hydrodynamic disturbances at the nozzle exit, thus closing the resonant loop. The screech 
phenomenon can produce sound pressure levels in excess of 160 dB in the vicinity of the jet 
(Raman 2 ). This phenomenon has been linked to structural damage in aircraft engine exhaust nozzles, 
premature wear in protective coatings, and increased engine noise. Details of the screech process and 
recent advances are described in Raman 3 . In this paper our focus is restricted to the demonstration of 
a robust interferometer that can provide accurate data even in the presence of intense acoustic fields. 
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An interferometer by definition requires two coherent beams to interfere, and environmental 
vibration and thermal disturbances affect the reference arm as well as the object arm making it 
difficult to isolate the information of interest. 

To extract quantitative results from interferometers phase-shifting interferometry 
(Greivenkamp and Bruning 5 ) was developed. It works on the principle that an interferogram, i.e. a 
sinusoidal function generated by the coherent interference of two optical waves, represents an 
equation with three unknowns: a DC bias, contrast, and phase. Therefore any of these three 
unknowns can be solved provided at least three equations exist. Multiple equations are obtained by 
recording interferograms with known phase differences added to the phase between the two optical 
waves. These equations are represented by: 

i n ( x > y)= ! avg ( x > y) + ft( x » y)c°s[4» (x, y) + e n ] 0) 

where I„ is the intensity in the n th recorded interferogram, I. vg is the average intensity, ri is the pattern 
contrast, <j> is the phase difference between the two optical waves, 0 n is an added phase, and (x, y) 
indicates that the variable is defined at each pixel. The phase <J> is directly related to the gas density, 
and is therefore the quantity of interest. This phase can be calculated simply and robustly by setting n 
= 1, 2,3,4, 5 and 0„ = nn/2 (Hariharan 6 ); knowledge of I« vg and ti is then not required provided they 

remain constant over the sequence of interferograms. 

The challenge is to generate the required sequence of Jt/2 phase steps. This is straightforward 
for conventional two beam interferometers like the Mach-Zehnder; typically a mirror in one beam is 
mounted on a piezoelectric device to precisely change the optical pathlength. However, two-beam 
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interferometers require a reference optical path, which must not experience flow disturbances. In 
large facilities this reference path can be quite long, and in noisy environments it is difficult to isolate 
even short paths from air turbulence. Single interferograms of jets have been recorded using short 
exposures and small facilities (Nakamura and Iwamoto 7 ) but fringe-tracking or Fourier transform 
techniques must be used for data processing. Common-path interferometers (Linnik 4 ) have been used 
to reduce the effect of environmental noise; they work because the reference beam travels the same 
path as the object beam so noise affects both beams in the same way. Since interferometry works by 
detecting phase differences between two beams, common disturbances do not affect the result. A 
modified point diffraction interferometer, the Liquid Crystal Point Diffraction Interferometer 
(LCPDI) was developed to both generate phase steps for automated data processing and provide the 
simplicity and robust performance of a common-path interferometer (see Fig. 2). This interferometer 
was developed for space flight experiments (Mercer and Creath 8 ) and is applied here to the case of a 
supersonic shock-containing jet. The LCPDI and the experiments are described, a density 
distribution measured simultaneously using the LCPDI and a Mach-Zehnder is presented, and flow 
measurements from LCPDI data is presented. 

H. EXPERIMENT 

The integrated line-of-sight density of a free jet was measured simultaneously with the LCPDI 
and a Mach-Zehnder interferometer. The experimental apparatus is shown in Figure 3. The optical 
path unique to the Mach-Zehnder is shown as double-hatched lines, the optical path unique to the 
LCPDI is shown single-hatched, and the path common to both interferometers is shown crosshatched. 
A frequency-doubled diode-pumped YAG laser operating at 532 nm was attenuated to 5mW by a 
neutral density filter (ND1) and directed by mirror Ml to a beam expander (EXPND) which generated 
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a collimated beam. Beamsplitter BS1 transmitted 50% of the beam to form the Mach-Zehnder's 
reference arm. The reflected beam was directed by mirror M3 to pass just downstream of a 4.5 mm 
ID circular nozzle. Beamsplitter BS2 transmitted 50% of this beam to form the Mach-Zehnder's 
object beam, which was recombined, with the reference beam by beamsplitter BS3. These beams 
were recorded by camera CAM2. The beam reflected by BS2 was focused by a f72.6 Cooke Triplet 
lens nearly onto the LCPDI. The LCPDI consists of two glass plates with a 9-micron gap filled with 
nematic liquid crystals and 9-micron microspheres. The LCPDI is positioned so that the focused 
beam illuminates one microsphere. (In practice several microspheres may be present in the focused 
beam.) Some of the light passes around the microsphere through the liquid crystals, and some of the 
light is diffracted by the microsphere forming a spherical wave. Both the diffracted reference wave 
and the "pass-through" object wave travel co-linearly to camera CAM1 where they combine 
coherently to form an interferogram. 

Phase stepping is accomplished in the Mach-Zehnder by applying a sequence of DC voltages 
to a piezoelectric element in mirror M2. This changes the path length of the reference beam (hence 
its phase) while leaving the object beam unaffected. Phase stepping is accomplished in the LCPDI by 
applying a sequence of AC voltages across indium tin oxide electrodes deposited on the inner 
surfaces of the glass plates. This changes the refractive index of the birefringent liquid crystals and 
steps the phase of the light passing through them, while leaving the phase of the diffracted beam 

unchanged. 


In both cases, the phase 4> is measured by recording a sequence of 5 interferograms with an 
additional id! radian phase step added each time. Each interferogram was recorded with an exposure 
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time of 1/30 sec and a delay of 2 seconds between frames to ensure that the liquid crystals were 
properly oriented. The data acquisition time was therefore 8.3 sec. A sequence of 5 frames was 
recorded with the jet on and then again with the jet off to provide a background measurement. A 
phase map was calculated from each of the 2 sets of 5 frames, and the two-phase distributions were 
then subtracted to produce a single array containing the optical phase difference between the flow and 
no-flow conditions. This array had the same number of elements as the recording camera, 752 x 480 
pixels in our case. Our field of view was 20 x 12.7 mm for a spatial resolution of 0.026 mm/pixel. 

The gas density, p, is determined from this measured phase difference, A<|> by: 

J p(x,y)ds = -^— A<j>(x,y) + Lp 0 ( 2 ) 

0 27 uv 

where X is the laser wavelength, K is the Gladstone-Dale constant, L is the physical path length 
through the jet, and p 0 is the no-flow density. For the data presented here, X=532.0 nm, K=0.2270 
cm 3 /g (Vest 9 ), L=50 cm and po =1.205 kg/m 3 . 

If the flow under study is two-dimensional, then the above integral simplifies to Lp and the 
density is simply calculated. If the flow is axisymmetric, then an Abel transform can be used to 
calculate the radial distribution of density at each axial location. The calculation of this transform 
using data from the LCPDI is greatly simplified because the high spatial density of the data 
eliminates the need for interpolation (Kassemi et. al. 10 ). If a folly three-dimensional flow is of 
interest, then data must be taken from multiple views and reconstructed tomographically. 
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III. EXPERIMENTAL UNCERTAINTY 

The uncertainty of the measured density is determined as follows. Since the density itself is 
extracted in one of the three ways outlined above, only the uncertainty in the line-of-sight integrated 
density is calculated. Define the relative, integrated density Q as: 


Q = J pds - Lp 0 ( 3 ) 

The (x,y) notation has been dropped for convenience. It is readily seen from equation 2 that Q is 
proportional to X and A*, and inversely proportional to K. The relative uncertainty in Q is therefore 

equal to: 



+ 




( 4 ) 


For ax=0.1 nm, a K =0.0001 cm 3 /g and the values for X and K given above, the first two squared terms 
equal 0.02% and 0.04%, respectively. Because these uncertainties in the wavelength and Gladstone- 
Dale constant are so small, the density uncertainty approximately equals the phase measurement 
uncertainty for phase measurement uncertainties greater than 1 degree. Therefore, if the phase is 
measured to 1 degree, 0^=1/360 and oq/Q=0.3%. The relative density uncertainty is <1% for 
measured phase uncertainties <3.6 degrees, and the density is known to <10% if the measured phase 
is determined within 36 degrees. The typical standard deviation for the measured phase is 4 degrees 
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and the data presented here spanned about 6 interference cycles. Uncorrelated speckle from frame-to- 
frame and other optical interference can cause noisy interferograms which increase the measured 
phase uncertainty. Note that this analysis holds for any single-pass interferometer. 

IV. RESULTS 

Figure 4 shows a spark-schlieren photograph of the intense acoustic environment produced by 
an underexpanded jet at a fully expanded Mach number (Mj) of 1.4. A larger scale rectangular nozzle 
with exit dimensions of 8.89 cm x 0.635 cm was used for this visualization. The view shown in 
Figure 4 is that of the narrow dimension of the nozzle. The oscillatory motion of the entire jet and 
feedback shocks in the nearfield is clearly visible. Sound pressure levels in the nearfield for such jets 
are generally in excess of 160 dB in the nearfield. Note also that the feedback shock waves 
(occurring outside the flow region) appear to originate at the end of the third shock-cell. 

In the rest of this section we will focus on the advantages of using the LCPDI technique to 
study such flows that produce a highly disturbed environment. 

For the interferometry measurements two convergent nozzles were used: a circular nozzle with 
an inner diameter of 4.5 mm, and a rectangular nozzle with inner dimensions 1.27 x 6.10 mm. The 
rectangular nozzle was oriented with its long axis parallel to the optic axis. Compressed air was used 
as the gas; its total pressure was 32.0 psig (ambient pressure was 14.5 psig) for the circular nozzle 
and 60.8 psig for the rectangular nozzle. There was a substantial pressure loss in each nozzle, so 
these pressures corresponded to 23.56 psig and 53 psig at the exit plane for the two jets as determined 
by a later measurement of exit pressure vs. tank pressure using a pitot probe (0.8 mm o.d.) at the exit 
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plane. The fully expanded jet Mach number was 1.26 and 1.66 for the circular and rectangular 
nozzles, respectively. 


Figure 5 shows the results for both the LCPDI and the Mach-Zehnder interferometers. The 
gray scale images are coded so that darker grays indicate lower density. In both cases three shock- 
cells are clearly present. The LCPDI has the uniform background that one would expect for this flow. 
The background in the Mach-Zehnder data shows a periodic structure that is charactenstic of phase 
stepping errors indicating that the reference beam did not remain unchanged during the data 
acquisition. This figure clearly shows the improved performance of the LCPDI's common path 
design. The data reduction was automated by the phase stepping process. 

The integrated density distribution represented by Eqn. 2 for the circular and rectangular 
nozzles are shown in Figure 6 (a) and (b). No smoothing has been done on this data, but the images 
were cropped down to 500x350 pixels to more clearly show the features of interest. Both images are 
shown with the same gray scale; the rectangular jet has a much larger density variation than the 
circular. The Mach diamonds are clearly visible in both images as the regions of darkest color (lowest 

density). 


Cross sections along the centerline from both data sets are shown in Figure 7. This data was 
averaged over a neighborhood of 7 pixels. Again, this is the line-of-sight integral of density divided 
by optical path so the vertical axis is best interpreted as an arbitrary scale. Both data sets are plotted 
to the same scale. The flat line at the left in both cases represents the nozzle with flow from left to 
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right. From this plot the shock spacing for each nozzle can be determined. The spacing for the 
rectangular nozzle is 4.40 mm and for the circular nozzle 4.32 mm. 

Seiner & Norum" measured the shock-cell spacing from circular nozzles that had various 
design mach numbers and provided a simple correlation for their data. Their correlation agreed 
reasonably with the theoretical wavelength derived by Pack . For a circular jet with a design Mach 
number of 1, Norum & Seiner 11 showed that the average shock spacing can be represented by L/D = a 
(Mi 2 ' 1) M > where Mj is the fully expanded jet Mach number, a — 1.16 and b = 1.0. The shock-cell 
spacing (L/h) for rectangular jets can be calculated using models developed by Tam 13 and Morris et 
al. 14 . Tam 13 used the Prandtl approach (Pack 12 ) and considered the shock-cell system in the jet 
column to be bounded by a mixing layer that was approximated by a vortex sheet. The problem was 
then solved by eigenfunction expansion and the solution to the eigenvalue problem yielded explicit 
shock-cell spacing. Morris et al. 14 also used the vortex sheet model but used a boundary element 
capable of accommodating nozzles with arbitrary cross-sections to obtain the results. A comparison 
between the present measured values and those of Tam 13 and Morris et al. 14 is given in Table I. Note 
that such linear shock-cell models are only expected to be accurate at low levels of underexpansion 
before the end of the first shock-cell turns into a Mach disk. 

V. CONCLUDING REMARKS 

A new, robust, phase-stepping, common-path interferometer was used to measure the density 
distribution in a free jet. This interferometer has superior vibration insensitivity compared to the 
Mach-Zehnder interferometer because of its common path design. This is useful for vibratory 
environments where a reference beam has to travel long distances and can be easily disturbed. The 
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present technique uses a microsphere in a liquid crystal layer to locally generate a reference beam that 
is not easily disturbed. The LCPDI is designed to work with a single microsphere, but in practice that 
has been difficult to achieve. Instead, several microspheres are present, causing extraneous 
interference. Also, dye is required in the liquid crystal layer to attenuate the object beam. Currently 
this dye changes its attenuation with applied voltage, and requires more complex data processing than 
standard phase stepping interferometers. Work is underway to correct these defects, but even with 
these drawbacks the data presented here shows that the LCPDI produces useful results. 

Although previous interferometry techniques had the potential to produce quantitative density 
data they were not very widely used because fringe-tracking methods are often tedious, subjective, 
and prone to errors. The present work uses a LCPDI with the phase shifting technique to generate 
information required to solve three different equations for the unknown phase distribution. In 
addition, automatic data acquisition and data reduction eliminate the "trained eye" requirements of 

fringe-tracking methods. 

Underexpanded supersonic jets were selected for the present study because they produce very 
intense acoustic tones with near-field levels in excess of 160 dB. In addition, the underexpanded jet 
also displays very organized and drastic density changes that were resolved in detail by the present 
technique for both rectangular and circular geometries. Shock spacing were also measured for two 
nozzles and the results compared favorably to theory. 

The LCPDI is thus an easy to use interferometer that works in harsh environments, is easy to 
align and is capable of measuring optical wavefronts with high accuracy. This technique has the 
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potential for very wide use for density measurements in research and industrial situations. Because 
the data acquisition time was nominally 8 sec, only steady-state flows can currently be measured. It is 
believed that this acquisition time can be reduced to 0.5 sec/frame by judicious application of 

stepping voltages. 
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Nozzle 

Total Pressure 
(psig) 

Normalized Shock Spacing 
(L/D, L/h) 

% Difference 

Circular 

Rectangular 

Rectangular 

23.56 

53.0 

53.0 

Experimental 

Calculated 

6.25 

0.43 

3.90 

0.96 ± 0.02 
3.46 ±0.06 
3.46 ±0.06 

0.90 (Seiner) 
3.475 (Tam) 
3.325 (Morris) 


Table I. Calculated and Measured Shock-Spacings. 




CAPTIONS FOR FIGURES 

Figure (1) Typical narrowband spectrum for an underexpanded jet showing the three distinct 
components of supersonic jet noise. 

Figure (2) Schematic of Liquid Crystal Point Diffraction Interferometer (LCPDI). 

Figure (3) Experimental configuration for simultaneous LCPDI and Mach-Zehnder data 
acquisition. 

Figure (4) Spark-schlieren photograph showing an underexpanded screeching jet at Mj — 1.4. 

Figure (5) Line-of-sight integrated density distribution downstream of a circular jet 

measured with (a) LCPDI (b) Mach-Zehnder interferometer. 

Figure (6) Measured line-of-sight integrated density for (a) circular and (b) rectangular jets. 
Flow is from left to right. 

Figure (7) Line-of-sight integrated density along center plane for circular and rectangular 
underexpanded jets. 
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Figure 1 . — Typical narrowband spectrum for an underexpanded jet showing the 
three distinct components of supersonic jet noise. 
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Figure 2. — Schematic of Liquid Crystal Point Diffraction Interferometer (LCPDI). 


E-1 1 356 Raman 9pt/1 00% J J 




Figure 3. — Experimental configuration for simultaneous LCPDI and Mach-Zehnder data 
acquisition. 
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Figure 4. — Spark schiieren photograph showing an underexpanded screeching 
jet at Mj =5 1 .4. 
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Figure 5.— Une-of-slght integrated density distribution downstream of a 
circular Jet measured with (a) LCPDI. (b) Mach-Zehnder interferometer. 





Figure 6.— Measured line-of-sight integrated density for (a) circular and (b) rectangular jets. 
Flow is from left to right. 
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Figure 7. — Line-of-sight integrated density along center plane for circular and rectangular 
underexpanded jets. 
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